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SUMMARY 

The  method  of  characteristics  is  used  to  present  general  compatibility  relations  for 
nonsteady  liquid  flow  or  water-hammer  theory  which  permit  the  investigation  of  the  dy¬ 
namics  of  the  flow  imder  diverse  conditions;  that  is,  for  a  liquid  flowing  in  a  variable 
area  conduit  affected  by  heat  addition,  by  gravity,  by  friction  forces,  and  by  centrifugal 
forces  such  as  in  a  mixed- flow  radial  vane  piimp.  The  compatibility  equations  are  pre¬ 
sented  in  finite  difference  form,  and  numerical- graphical  methods  for  their  solution  are 
indicated. 


INTRODUCTION 

In  rocket  engine  and  Rankine  cycle  systems ,  the  analysis  of  nonsteady  liquid  flow  is 
often  complicated  by  the  transmission  of  the  fluid  in  a  conduit  with  a  gradually  varying 
cross-sectional  area.  In  addition,  the  fluid  may  be  affected  by  heat  addition  and  a  num¬ 
ber  of  forces  such  as  gravity,  dissipative,  and  centrifugal  (e.  g. ,  as  in  a  pump).  Hence, 
a  general  method  is  required  which  will  allow  the  investigation  of  fluid  transients  under 
such  diverse  conditions. 

Analyses  for  nonsteady  liquid  flow  have  been  presented  in  a  niunber  of  papers  for 
particular  conditions;  that  is,  for  a  tapered  line  (ref.  1),  for  friction  (ref.  2),  and  for  a 
pump  as  a  lumped  component  in  a  conduit  (ref.  3). 

In  this  report  the  theory  of  nonsteady  liquid  flow  as  developed  in  reference  4  will  be 
extended  so  that  the  dynamics  of  a  liquid  in  a  conduit  with  a  gradually  varying  cross- 
sectional  area,  undergoing  heat  addition,  and  affected  by  the  previously  mentioned  forces 
may  be  analyzed.  Applications  will  be  outlined  for  nonsteady  liquid  flow  with  heat  addi¬ 
tion,  for  tapered  conduits  with  gravity  and  dissipative  forces,  and  for  the  internal  dy¬ 
namics  of  a  mixed  flow  pump  with  radial  vanes. 


ANALYSIS 


Derivation  of  Fundamental  Compatibility  Relations 


Continuity.  -  The  continuity  equation  (ref.  5)  is  given  by 


—  (pA)  +  —  (puA)  =  0  (1) 

at  a? 


O  A 

and  if  —  =  0,  equation  (1)  may  be  expanded  to 

at 

Aiet/ieuH.i!!p\AK-“pu  =  o 
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The  term  — ,  a  function  of  the  elasticity  of  the  conduit,  is  omitted  since  it  is  usually 
^t  . 

negligibly  small  compared  to  the  variation  of  area  with  distance  along  the  conduit  — . 

H  equation  (2)  is  divided  by  pA,  then 


1  +  i  u  ^ +  u  =  0 

p  at  p  a|  a|  a| 


and  since 


(3a) 


P  =  f(p,  s)  (3b) 

equation  (3)  may  be  expressed  by 

pVap  at  as  at/  pVapai  as  a^/  a|  a? 

Assuming  that  is  negligibly  small  for  a  liquid  and  using  the  definition  for  the  acoustic 

as 

velocity,  namely. 


(5) 


result  in 
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(6) 


1  3P,  1  9P„^au^^^3(lnA)_n 

pa2  pa2 

Inasmuch  as  the  pressure  P  is  expressed  in  force  per  unit  area  and  the  density  p  in 
mass  per  imit  volume  and  the  head  in  the  units  of  pressure  divided  by  the  specific  weight 
of  the  fluid  in  the  Earth’s  gravitational  field,  the  following  relation  may  be  deduced 

P  -  P 

- ^  =  geH  (7) 

where  P  is  the  datum  pressure,  a  constant.  In  view  of  equation  (7)  and  the  assumption 
that  the  variation  of  the  density  has  a  negligibly  small  effect  on  the  head  (i.  e. ,  p  =  const. ), 
equation  (6)  may  be  presented  by 

gc  8H  ,  gc^8H^3u^^^3(lnA)_o  (8) 

^2  at  ^2  a|  a|  a| 

Multiplying  by  a  /g  and  transposing  some  terms  result  in 

an  ,  ^  aH  ,  a^  au  _  ua^  a  (in  A)  (9^ 

at  ai  gg  ai  g^.  a| 

Momentum.  -  Consider  the  momentum  equation  given  by  (ref.  5) 

iH  +  uiH  =  -l^  +  f  (10) 

at  a|  p  a| 

Multiplying  equation  (10)  by  a/g^  and  considering  equation  (7)  give 

l_ilH+i_u^  +  aiS=i_f  (11) 

gc9t  gc  9^  91  g^ 


Subtracting  equation  (11)  from  equation  (9)  yields 

i- (h  .  A  uU  (u  .  a)  iS  .  (u  -  a  iH  \  . 

at  \  g  /  a^  g  \  a|  a^/ 


ua^  a  (In  A)  a  ^ 

g„  a?  g„ 
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if  S’/ Sf,  is  a  constant.  By  further  algebraic  manipulation,  equation  (12)  can  be  pre¬ 
sented  in  the  form 


(n-iu) +(u-a)±/H-l-u'\  -  -  ■*)  »  f 

V  Sc  /  V  J  g^  H  g. 


Adding  equations  (9)  and  (11)  yields  another  expression,  namely. 


(h  uK  (i,  ..  a)  A  fw  ..  A  iA  -  -  9(ln  A)  ^  a  ^ 

St  \  gc  /  9  A  gp  /  g,  H  g„ 


Now  let 


H  +.^u  =  C 


H  -  u  =  C 


and  define  two  directional  derivatives  as  follows: 


6t  at  ai 


£=i.*(u.a)A 

6t  at  a^ 


With  the  help  of  equations  (15)  to  (18),  equations  (13)  and  (14)  may  be  expressed  by 


ua^  a(ln  A)  a 

gc  9^  gr 


_  _  ua^  a(ln  A)  ^  a  ^ 


Sc 
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In  the  terminology  of  the  method  of  characteristics  (ref.  5,  vol.  1,  appendix  A),  equa¬ 
tions  (19)  and  (20)  can  be  considered  as  right  running  and  left  ruiming  compatibility  equa¬ 
tions  respectively,  and  their  characteristic  directions  in  the  time-distance  plane  are 


«=u-a 


along  the  right  running  characteristic  and 

^  =  u  +  a  (22; 

dt 

along  the  left  ruiming  characteristic.  The  similarity  of  these  results  and  those  of  refer¬ 
ence  4  is  evident. 

If  the  flow  velocity  is  very  small  compared  to  the  acoustic  velocity,  equations  (21) 
and  (22)  become  simply 

^  =  -a  (23 


If  the  two  terms  appearing  to  the  right  of  equations  (19)  and  (20)  are  omitted  or  set  equal 
to  zero,  then 


1  =  0 


(25a) 


C  =  H  -  —  u 

y  or 


(25b) 


is  the  right  running  compatibility  equation  and 


Cx  =  H+-u 
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Is  the  left  running  compatibility  relation.  These  equations  are  the  more  familiar  water- 
hammer  relations. 

Energy.  -  The  energy  equation  gives  rise  to  a  third  compatibility  relation  and  char¬ 
acteristic  direction.  This  may  be  seen  by  considering  the  energy  equation  in  a  conduit  as 
expressed  by  (ref.  5) 


qpA  =  —  pA/i  +  — )  +—  puA(i  +  Pv  +  — 

0t  \  2  /  0^  I  \  2 


which  after  some  manipulation  becomes 


V  2/L  at  05  J  05  Vat  05/  L0tV2/  05  V 2 


By  the  definition  of  the  co-moving  or  substantial  derivative,  namely, 


Dt  0t  dt  05 


(29a) 


D=i-  +  uA 

Dt  0t  05 


(29b) 


and,  in  view  of  equation  (1),  equation  (28)  may  be  presented 


p  +  uA  —  +  pA  —  +  pAu  —  =  qpA 
05  05  Dt  Dt 


From  continuity,  the  following  relation  is  obtained: 


0(uA) 


J^^  =  Ap_A£ 
05  Dt 


where  —  is  omitted.  Hence,  instead  of  equation  (30), 

at 
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Time, 


+  pA  —  =  pAq 


(32) 


it  follows  that 


(35) 


if  the  magnitude  of  the  term  uf,  which  represents  in  the  negative  sense  the  specific  work 


Distance  along  conduit,  ^ 


Figure  1.  -  General  characteristic  net  point  for  nonsteady 
liquid  flow  or  water  hammer. 


done  by  dissipative  forces  per  imit  time,  is  much 
smaller  than  the  value  of  q  (i.  e. ,  the  heat  addi¬ 
tion  is  reversible).  Equation  (35)  is  the  third 
compatibility  equation;  and  in  view  of  equations 
(29a)  and  (29b),  the  direction  of  the  characteris¬ 
tic  in  the  time -distance  plane  will  be  given  by 

^  =  u  (36) 

dt 

Hence,  in  general,  a  net  point  in  nonsteady  liq¬ 
uid  flow  is  defined  by  three  characteristics  j  the 
projection  of  which  on  the  time-distance  plane 
is  seen  in  figure  1. 

Heat  addition.  -  Equations  (15)  and  (16)  are 
two  equations  with  two  unknowns,  H  and  u, 
since  C  is  a  constant  along  the  left  running 
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characteristic  and  Cy  is  a  constant  along  the  right  running  characteristic,  and  the 
acoustic  velocity  a  is  also  constant.  If  none  of  the  parameters  of  equations  (15)  and  (16) 
are  affected  by  entropy  change,  which  can  be  reasonably  assumed  in  liquids,  the  contri¬ 
bution  of  heat  addition  to  the  dynamic  analysis  of  liquids  may  be  omitted.  Yet  it  is  desir¬ 
able  to  estimate  the  temperature  variation  of  a  liquid  as  it  proceeds,  for  example, 
through  a  boiler  or  heated  pump  inlet  so  as  to  predict  the  incipience  of  boiling.  This  may 
be  accomplished  by  following  a  particle  of  fixed  identity  in  the  fluid  as  implied  in  equa¬ 
tion  (35).  From  thermodynamics,  it  is  given  that  for  constant  volume  reversible  heating 


T  Ds  =  CyDT 

Integrating  equation  (37)  along  the  particle  path  (cf.  fig.  1)  results  in 


T3  =  T4e 


(S3-S4)/C^ 


and,  consequently,  by  equation  (35) 


(37) 


(38) 


^3  = 


’;'(t3-t4)/C^ 


where 


(39) 


T 


Evaluation  of  the  Specific  Force 

The  specific  force  f  (or  force  per  unit  mass)  appearing  in  the  fundamental  momen¬ 
tum  equation  (eq.  (10))  is  undefined.  It  may  be  equated  to  a  specific  gravity  force,  or  a 
dissipative  force,  or  a  specific  centrifugal  force,  or  the  sum  of  all  these  depending  on 
the  particular  study  to  be  made.  In  the  following,  the  specific  force  f  will  be  equated, 
first,  to  the  sum  of  the  specific  gravity  and  specific  dissipative  forces  and,  then,  in  order 
to  analyze  the  internal  dynamics  of  a  mixed-flow  radial  pump,  to  the  sum  of  specific  cen¬ 
trifugal  and  specific  dissipative  forces. 

Gravity  and  dissipative  forces.  -  If  the  gravity  force  per  unit  mass,  or  in  other 
words  the  specific  gravity  force  g,  is  opposite  in  direction  to  the  assumed  positive  direc¬ 
tion  of  the  specific  force  f,  the  sum  of  the  specific  gravity  and  dissipative  forces  may  be 
expressed  by 
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(40) 


(46) 


where  the  subscripts  refer  to  net  points  in  the  characteristic  network  (cf.  fig.  1).  The 
subscripts  13  and  23  represent  average  values  between  net  point  1  and  3  and  net 
point  2  and  3,  respectively.  Clearly,  equations  (43)  and  (44)  may  be  more  easily  eval¬ 
uated  if  the  variation  of  area  with  conduit  length  is  expressed  in  the  form 
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A  =  e' 


(47a) 


e(ln  A)  ^  Sh(|) 
31  31 


(47b) 


Figure  2.  -  Mixed-flow  pump  with  radial  vanes.  Seeing  that 
F  ■  5  sin  0,  it  follows  that  f  =  sin^  0  +  ij?rj  sin  0. 


vanes  may  be  expressed  by  (cf.  fig.  2) 


where  h(^ )  is  some  easily  differentiable 
fimction  of  conduit  length. 

Pump  dynamics.  -  The  equation  of 
motion  for  a  mixed-flow  pump  with  radial 


9u  9u  /  2  V  .  2  2  \ 

—  +  u - (o)  I  sm  a  +  w  r.  sin  a)  = 

9t  0  f  ^  / 


P  9^ 


where  u  is  the  through-flow  component  of  the  flow  velocity  relative  to  the  blade  and  a 
is^the  angle  the  through-flow  component  makes  with  the  axis  of  the  pump.  The  terms 
cj  ^  sin  CT  -I-  ct)  r.  sin  o  represent  the  centrifugal  force  per  unit  mass.  H,  in  addition, 
friction  forces  are  significant,  then 


f  =  sin^CT  -I-  w^r.  sin  a  -  — H 

^  2D  Iu 


and  equations  (41)  and  (42)  may  be  modified  to 


gp  3|  g 


y_  ua'^9(lnA)  a.  (..2 2  2  n 

- ^ - - - w  Ism  a  +  w  r.  sm  a  -  — — 


2D  u 


«  3|  g, 


_  ua*^  9(lnA)  ,  a  /  2^  .  2  2  .  u 

- ^ - -  +  —  w  I  sm  a  +  w  r.  sin  a  -  --iL  _!i 


2D  u 


or  in  finite  difference  form  to 


‘^>'3  '  ^  “3  -  ^  a,2rj  sto  ^  (tj  -  tj) 

23  L  /_  00 
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#3  -  tj) 


r 


ua^  0(ln  A) 


L°C 


(t3  -  tj)  + 


13 


2^  _.„2 


cl  41  V  •  ^  CJ 

_  a>  I  sin  a  +  CO 


cr 


Sin 


Xu^  u 
2D  |u 


13 


(52) 


If  cavitation  occurs  in  the  pump  (the  head  is  less  than  the  vapor  head),  a  vapor  region 
will  appear  in  the  characteristic  network.  The  vapor  head  and  the  characteristic 
will  define  the  flow  into  the  cavitation  zone,  and  the  Cy  characteristic  will  define  the 
flow  out.  The  time  increment  and  past  conditions  will  help  to  estimate  the  augmentation 
in  the  size  of  the  cavitation  zone.  Further  details  on  the  analysis  of  the  dynamics  of  cavi¬ 
tation  regions  by  essentially  the  method  of  characteristics  may  be  found  in  reference  6. 


CONCLUDING  REMARKS 

General  compatibility  relations  have  been  presented  which  will  permit  the  analysis  of 
the  dynamics  of  a  liquid  under  diverse  conditions.  In  view  of  their  complexity,  a  high¬ 
speed  computer  should  be  employed  in  their  practical  application  to  any  particular 
problem . 

Although  an  electromagnetic  pump  has  not  been  considered  in  this  memorandum,  it 
is  clear  that  in  place  of  the  centrifugal  force  per  unit  mass  term  in  the  compatibility  equa¬ 
tions  appertaining  to  pump  dynamics,  an  electromagnetic  force  per  imit  mass  term  may 
be  employed. 

The  characteristic  methods  of  this  report  may  be  applied,  for  example,  to  an  approx¬ 
imate  estimate  of  the  effect  of  the  through-flow  channel  of  a  pump  on  flow  dynamics .  Fur¬ 
ther,  the  techniques  of  this  report  may  be  used  to  predict  the  dynamic  formation  and  col¬ 
lapse  of  cavitation  regions  in  a  heated  pump  inlet  line  as  well  as  the  incipience  of  two- 
phase  flow  in  a  boiler  under  dynamic  conditions. 

Lewis  Research  Center, 

National  Aeronautics  and  Space  Administration, 

Cleveland,  Ohio,  Jime  8,  1966, 

120-27-04-27-22. 
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APPENDIX  A 


SYMBOLS 


A  area,  sq  ft 

a  acoustic  velocity,  ft/sec 

Cy  specific  heat  at  constant  vol¬ 

ume,  ft-lb/(slug)(°R) 

Cjj  parameter  (see  eq.  (15)),  ft 

Cy  parameter  (see  eq.  (16)),  ft 

D  hydraulic  diameter,  ft 

f  specific  force  or  force  per 

2 

imit  mass,  ft/sec 

g  specific  gravity  force,  ft/sec^ 

gp  acceleration  due  to  Earth's 

gravitational  field,  ft/sec 

H  head,  ft 

i  internal  energy,  ft-lb/slug 

P  pressure,  Ib/sq  ft 

q  heat  added  per  unit  mass  per 

unit  time,  ft-lb/(slug)(sec) 

r  radial  height  (see  fig.  2) 

r^  inlet  radial  height  (see  fig.  2) 


s  specific  entropy,  ft- lb/ 

(slug)(°R) 

T  temperature,  °R 

t  time,  sec 

u  flow  velocity,  ft/sec 

V  specific  volume,  cu  ft/slug 

A.  friction  coefficient 

^  distance  along  conduit,  ft 

p  density,  slug/cu  ft 

CT  angle  that  the  relative 

through-flow  component  of 
velocity  makes  with  pump 
axis 

parameter,  q/T 

00  angular  speed,  sec"^ 

Subscripts: 
o  datum 

1, 2,  3, 4  refers  to  net  point  locations 
(see  fig.  1) 
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